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ABSTRACT 

In order to investigate the origin of quasars, we estimate the bias factor for low-luminosity quasars 
at high redshift for the first time. In this study, we use the two-point angular cross-correlation function 
(CCF) for both low-luminosity quasars at —24 < M1450 < —22 and Lyman-break galaxies (LBGs). 

Our sample consists of both 25 low-luminosity quasars (16 objects are spectroscopically confirmed low- 
luminosity quasars) in the redshift range 3.1 < z < 4.5 and 835 color-selected LBGs with < 25.0 
at 2: ^ 4 in the COSMOS field. We have made our analysis for the following two quasar samples; 

(1) the spectroscopic sample (the 16 quasars confirmed by spectroscopy), and (2) the total sample 
(the 25 quasars including 9 quasars with photometric redshifts). The bias factor for low-luminosity 
quasars at 2 ~ 4 is derived by utilizing the quasar-LBG CCF and the LBG auto-correlation function. 

We then obtain the 86% upper limits of the bias factors for low-luminosity quasars, that are 5.63 and 
10.50 for the total and the spectroscopic samples, respectively. These bias factors correspond to the 
typical dark matter halo masses, log (MDM/(^~^AfQ)) =12.7 and 13.5, respectively. This result is 
not inconsistent with the predicted bias for quasars which is estimated by the major merger models. 
Subject headings: cosmology: large-scale structure of Universe — galaxies: active — quasars: general 
— surveys 
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1. INTRODUCTION 

The observed close relationship between the mass 
of the spheroidal component of a galaxy and its cen¬ 
tral supermassive black hole (SMBH) suggests that 
the evoluti on of galaxies and SMBHs are closely re¬ 
lated (e.g., IMarconi fc Himtl I2QQ3I: iHaring fc Rixl I2Q04I: 
iMcConnell fc Mai I2013D . Accordingly some important 
questions arise as when and how such a co-evolution has 
been established and what physical processes are essen¬ 
tially important. A straightforward approach to explore 
these issues is investigating the statistical properties of 
quasars at high redshifts, where the quasar activity (that 
corresponds to the growth of SMB Hs) is much more ac¬ 
tive than in the loca l universe fe.g.. lRichards et 311120061 
iGroom et al.ll2009al) . because different evolutionary sce¬ 
narios for SMBHs predict different statist ical properties 
of qu asars as function of redshift (e.g., I Hopkins et al.l 
[2003). 

The quasar activity is thought to be powered by mass 
accretion o nto a SMBH at the center of massive galaxies 
(iReedllO^ . The most efficient gas fueling mechanism is 
thought to be majo r and minor merg ers of gal axies (e.g., 
Sanders et al.lII988I iTaniguclil 119991 see also iTaniguclil 
201,111 . Therefore, in order to constrain the triggering 
mechanism for quasar activity, we need detailed studies 
of environmental properties of quasars. 

Motivated by these considerations, the two-point 
auto-correlation function (ACE) of quasars has been 
studied based on wide-field survey data, e.g., 2dF 
Quasar Red shift Survey a nd Sl oa n Digital Sky Sur¬ 
vey f e.g.. iCroom et ni 120051: iPorciani &: Norbe'r^ 
IMO^ iM^rs et al.l I20(M I2007t IShen et al.l IIoOtF 
Ida Angela et al.l 120081: IShen et HI 120091 : IRoss et ^ 
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I2009t llvashchenko et al] l2010t iWhite et alJ I2012D . The 
ACFs have also been investigated for active galactic 


Cilli et al.l 120091: iStarikova et al. 

120111: iKrnm 

oe et al.l 

20121: iKoutoulidis et al.l 120131: 

Allevato et al. 

12 QM). 


Most of these studies have shown that luminous 
AGNs tend to live in massive dark matter halos 
(~ 10^^ — suggesting that luminous 

AGN activity is triggered by galaxy mergers because 
such massive dark matter halos could be assembled by 
successive mergers of small dark matter halos. It is 
also reported that galaxy mergers do not account for 
the maj ority of the moderat e X-ray luminous AGNs at 
z < 2.2 (|Alievato et al.l[201111 . However, it is difficult to 
investigate small-scale clustering properties of quasars 
because of their low number density. 

Here it should be noted that the mass of dark-matter 
halos hosting quasars can be estimated also through 
the two-point cross-correlation function (CCF) for 
quasars and galaxies, combined with the AGF for 
galaxies. The advantage of the CCF is that the required 
size of the quasar sample is relatively smaller than 
the AGF analysis, though we need enough number 
of galaxies around quasars for the CCF analysis. In 
addition, it is possible to study the small-scale clustering 
properties of quasars through the CCF, that is too 
challenging for the ACF study. Therefore, in order to 
understand the triggering mechanism of quasar activity, 
it is useful to investigate the CCF for quasars and 
galaxies around them. Sev eral pioneering studie s have 


been made to date (e.g.. lAdelber^er & Steidell 

20051 

Padmanabhan et al.l 120091: Mountrichas et al.l 

2 ^ 

Coil et al.l 20091 iMiyaii et al.l 20111 Hickox 

et al. 

20111 iZhang et al.l 20131 iShirasaki et al.l 

2011 ; 

Mountrichas & Ceorgakakid 120121 IKomiva et al.l 

2013; 

Shen et al.1 120131 iKrolewski & EisensteinI 1201511. 
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instance, Zhang et al. (2013) investigated the spatial 
clustering of galaxies around quasars at redshifts from 
0.6 to 1.2. They found that the clustering amplitude 
is significantly larger for quasars with more massive 
black holes, or with bluer colors, while there is no 
dependence on quasar luminosity. This suggests that 
the mass of dark matter halos in which quasars reside is 
not correlated with the quasar luminosity. In addition, 
it is possible that the triggering mechanism of high- 
and low-luminosity quasars may be the same in this 
luminosity range. 

The CCF of quasars and galaxies has also been in¬ 
vestigated at higher redshifts. Shirasaki et al. (2011) 
investigated the projected CCF of AGNs and galaxies 
at redshifts from 0.3 to 3.0 and found significant excess 
of galaxies around the AGNs. They found that AGNs 
at higher redshifts reside in denser environments than 
those at lower redshifts. This suggests that major merg¬ 
ers are the preferred mechanism to trigger AGN activity 
at high redshifts. They also reported that there is no 
luminosity dependen ce of AGN clustering. At z ^ 3, 
iFrancke et al.l (j2008l l studied the two-point angular CCF 
of AGNs and Ly man break galaxies (LBGs) (see also 
iBielbv et al.ll201ll) and found that AGNs tend to be clus¬ 
tered more strongly than LBGs. They also found no lu¬ 
minosity dependence of AGN clustering. At z > 3, clus¬ 
tering of LBGs around quasars has been recently studied 


(e.g.. iHusband et ^1201311 . However the luminosity de¬ 
pendence of quasar clustering is not studied, due to the 
lack of adequate samples of low-luminosity quasars and 
galaxies around them. Therefore, the triggering mecha¬ 
nism of low-luminosity quasars at z > 3 has not yet been 
studied so far. 

Motivated by the situation described above, we study 
the two-point angular CCF of the 25 low-luminosity 
quasars (16 objects are spectroscopically confirmed low- 
luminosity quasars) and 835 LBGs in the redshift range 
3.1 < z < 4.5 in the COSMOS field. Then we derive the 
bias factor for low-luminosity quasars and constrain the 
dark matter halo mass in which low-luminosity quasars 
at z ^ 4 exist. The outline of this paper is as follows. In 
Section 2, we describe the data of low-luminosity quasars 
focused in this study and the method used for the photo¬ 
metric selection of LBG candidates. In Section 3, we re¬ 
port the results of the clustering analysis of quasars and 
LBGs. In Section 4 and 5, we give our discussion and 
summary. Throughout this paper we adopt a ACDM 
cosmology with 11^ = 0.3, Ha = 0.7, as = 0.9, and a 
Hubble constant of Hq = 70 km s“^ Mpc“^. We use the 
AB magnitude system. All of the errors reported in this 
paper are 1 sigma. 

2. THE SAMPLE 

2 .1. The Cosmic Evolution Survey 

Wide and deep multi-wavelength data are publicly 
available in the Cosmic Evolution Survey (COSMOS) 
field (Scoville et al. 2007). Therefore we can select 
large numbers of high-z galaxies and quasars in the 
same field to study their statistical properties such as 
the CCF. Another advantage of the COSMOS dataset 
is the dense sampling in the optical wavelength with 
intermediate-band Liters, that makes the estimates of the 
photometric redshifts of gala xies far more accu rate than 
in other deep-survey fi elds (lllbert et al.l 120091: see also 
iCardamone et al.l 12010(1 . For the above reasons, we de¬ 
cided to focus on the COSMOS field for the CCF analysis 
at z ^ 4. 

COSMOS is a treasury program of the Hubble Space 
Telescope (HST). It comprises 270 and 320 orbits al- 
located in the HS T Cycles 12 and 13, respectively 
(jScoville et al.ll2007t Koekemoer et al. 2007). The COS¬ 
MOS held covers an area of ^ 1.4° x 1.4° square which 
corresponds to ^ 2deg^, centered at R.A. (J2000) = 
10:00:28.6 and Dec. (J2000) =-L02:12:21.0. We use 
an upgraded version of t he ph otometric redshift cat¬ 
alogu e from lllbert etld] (1200911 (see also iCaoak et al.l 
l2007f) includi ng the new Ultr aVISTA data from the 
DRl release (iMcCracken et al.ll2012ll . to select samples 
of both quasars and LBGs at z ~ 4. This catalog 
covers an area of ~ 2 deg^ and contains several pho¬ 
tometric measurements. Specihcally in this paper, we 
use the u*-band 3 diameter aperture apparent magni- 
tude measured on t he image obtained with MegaCam 
(jBoulade et al.l[^003fl on the Canada-France-Hawaii Tele¬ 
scope (CFHT), and the 3 diameter a perture apparent 
magn itudes of the g'-, r'-, and z'-bands (|Taniguchi et al.1 
120071 1 measured o n the image obtained with the Sub¬ 
aru Suprime-Cam (jMivazaki et al.ir2002l l. The 5a limit¬ 
ing AB apparent magnitudes are u* = 26.5, g' = 26.5, 
r' = 26.6, and z' = 25.1 (3 diameter aperture). We 
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also use the Advanced Camera for Surveys (ACS) cata¬ 
log (jKoekemoer et al.ll2007l : iLeauthaud etalfc00711 when 
we select the low-lumin osity quasars to separate galaxies 
from point sources fsee llkeda et al.ir 201 lL I 20 T^ for more 
details). 

2 . 2 . Selection for Quasars at z ^ 4 

In order to calculate the bias factor for low-luminosity 
quas ars at z 4, we ne ed a low-luminosity quasar sam¬ 
ple. Ilkeda et al.l (|2011fl identified eight low-luminosity 
quasars through spectroscopic follow-up observations for 
their optical color-selected photometric quasar candi¬ 
dates. Additional low-luminosity quasars were also found 
by using the multi-wavelength imaging data includ¬ 
ing the optical broad, intermediate, narrowband, near- 
infrared and Spitzer /IRAC photometric measurements 
(|Masters et al.lf^012ll . To select low-luminosity quasars 
at z ~ 4 in the COSMOS field we used the following 
selection criteria: 

- 24.0 < Mi 45 o < -22.0, (1) 

and, 

3.1 <z <4.5. (2) 

Note that we adopt both z and M 1450 given in Masters et 
al. (2012). Then we reject objects which lie in masked 
regions (Capak et al. 2007). As a result, we selected 
25 quasars at 3.1 < 2 ; < 4.5. Among them, 16 objects 
are spectroscopically confirmed quasars. We refer this 
sample as the spectroscopic sample while the sample of 
25 quasars is as the total one. Total sample have been 
selected by utilizing the 29-band photometric data to re¬ 
move contaminants. Since they have been selected by 
utilizing such a lot of photometric data, we consider that 
the contamination rate for total sample is to be very 
low. Table 1, Figure 1, and Figure 2 show properties of 
these quasars, their redshift distribution, and their mag¬ 
nitude distribution, respectively. The median, mean, and 
the standard deviation of the redshift of the 25 quasars 
(16 spectroscopically confirmed quasar redshift) are 3.45 
(3.59), 3.59 (3.68), and 0.40 (0.40), respectively. 

2.3. Selection for Lyman Break Galaxies at z ~ 4 

We also select a sample of Lyman break galaxies 
(LBGs) in the redshift range 3.1 < z < 4.5 in the COS¬ 
MOS field, utilizing the two color diagram of r' — z' vs. 
g' — r' (Figure 3). For the selection of LBGs at z ^ 4, 
we adopt the following selection criteria: 



149.411400 < R.A.(degree) < 150.826934, 

(3) 


1.49878 < Decl.(degree) < 2.91276, 

(4) 


z' < 25.0, 

(5) 


u* > 27.05, 

(6) 


5 '-r'> 1 . 0 , 

(7) 


1—1 

+ 

I 

1—1 

A 

I 

(8) 

and. 

r' — z' < 1.5, 

(9) 

where u* = 27.05 corresponds to the 3cr limiting mag¬ 
nitude in the u*-band. The criterion (8) is adopted to 
select LBGs without significant contamination from low- 


z elliptical galaxies (see Figure 3). In order to remove 


TABLE 1 

Properties of the low-luminosity quasars at 2 ~ 4 


ID“ 

R.A. 

(deg) 

Decl. 

(deg) 

Afi450 

(mag) 

^sp 

z*’ 

adopt 

298002 

150.43706 

1.649305 

-22.34 

3.89 

3.89 

329051 

150.16891 

1.774590 

-22.71 

- 

4.35 

330806 

150.10738 

1.759201 

-22.85 

4.14 

4.14 

381470 

149.85396 

1.753672 

-22.11 

- 

3.30 

422327 

149.70151 

1.638375 

-22.49 

3.20 

3.20 

507779 

150.48563 

1.871927 

-23.78 

4.45 

4.45 

519634 

150.27715 

1.958373 

-22.61 

- 

3.40 

710344 

150.62828 

2.006204 

-22.06 

- 

3.45 

804307 

150.00438 

2.038898 

-23.56 

3.50 

3.50 

887716 

149.49590 

1.968019 

-22.38 

- 

3.23 

1046585 

149.85153 

2.276400 

-22.39 

3.37 

3.37 

1060679 

149.73622 

2.179933 

-22.23 

4.20 

4.20 

1110682 

149.50595 

2.185332 

-22.71 

- 

3.28 

1159815 

150.63844 

2.391350 

-22.98 

3.65 

3.65 

1163086 

150.70377 

2.370019 

-23.00 

3.75 

3.75 

1271385 

149.86966 

2.294046 

-23.42 

3.35 

3.35 

1273346 

149.77692 

2.444306 

-22.65 

4.16 

4.16 

1371806 

150.59184 

2.619375 

-22.10 

- 

3.12 

1465836 

150.13036 

2.466012 

-22.56 

3.86 

3.86 

1575750 

150.73715 

2.722578 

-22.78 

3.32 

3.32 

1605275 

150.62006 

2.671402 

-22.82 

3.14 

3.14 

1657280 

150.24078 

2.659058 

-22.74 

3.36 

3.36 

1719143 

149.75539 

2.738555 

-22.26 

3.52 

3.52 

1730531 

149.84322 

2.659095 

-22.15 

- 

3.51 

1743444 

149.66605 

2.740230 

-22.54 

- 

3.15 


“ ID for Table 3 of [Masters et all 11201211 . 

^ Redshift given by Masters et al. (2012) and adopted in our 
analysis. 

low-z objects, we add the criteria ( 6 ), (7), and (9). These 
selection criteria are adjusted to select LBGs with photo¬ 
metric redshifts (zph), whose distribution is similar to the 
redshift distribution of our quasar sample. We also re¬ 
move eight spectroscopically confirmed quasars satisfying 
the above criteria from the LBG sample, because those 
8 objects are not LBGs apparently. Note that seven ob¬ 
jects among the eight removed objects are included in our 
quasar sample while the remaining one object is not, be¬ 
cause its magnitude (M 1450 = —20.91) is out of the mag¬ 
nitude range of our quasar sample (—22 < M 1450 < —24). 

As a result, we obtain a sample of 835 LBGs at z ~ 4 
in the COSMOS field. We use the upgraded photomet¬ 
ric redshift (originally described in Ilbert et al. 2009 
with including the new UltraVISTA DRl data) to in¬ 
vestigate the photometric redshift distribution of the 
color-selected galaxies. Figure 4 shows the photometric 
redshift distribution of color-selected galaxies. The me¬ 
dian, mean, and the standard deviation of the redshift 
of the color-selected galaxies are 3.59, 3.13, and 1.32, re¬ 
spectively. As shown in Figure 4, there are some low-z 
(z < 1) galaxies in the color-selected LBG sample. We 
consider that most of these contaminants are low-z ellip¬ 
tical galaxies, being inferred from the color track of the 
model elliptical galaxy shown in Figure 3. Note that it 
is difficult to remove these low-z contaminants by mod¬ 
ifying the adopted color-selection criteria, because the 
optical colors of LBGs and those contaminants are sim¬ 
ilar. The median, mean, and the standard deviation of 
the photometric redshift of the color-selected galaxies, 
after removing objects whose photometric redshifts are 
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Fig. 1.— Redshift distribution of the low-luminosity quasars at 
3.1 < z < 4.5 used in this work. Upper and lower panels show the 
redshift distribution of 25 low-luminosity quasars (total) and 16 
spectroscopically confirmed low-luminosity quasars, respectively. 


Fig. 2.— Magnitude distribution of the low-luminosity quasars at 
3.1 < 2 : < 4.5 used in this work. Upper and lower panels show the 
magnitude distribution of 25 low-luminosity quasars (total) and 16 
spectroscopically confirmed low-luminosity quasars, respectively. 


below 3, are 3.67, 3.71, and 0.62, respectively. These re¬ 
sults are similar to that of low-luminosity quasars in this 
work. The typical error of the photometric redshift for 
the color-selected LBGs is ^ 0.1. There are 67 spectro¬ 
scopically confirmed objects in our color-selected galax¬ 
ies and their redshift distribution is shown in Figure 5. 
We confirm that the spectroscopic redshift distribution 
of color-selected galaxies is also similar with that of low- 
luminosity quasars in this paper. 

Figure 6 shows the magnitude distribution of the color- 
selected galaxies used in this paper and we also con¬ 
firm that the magnitude distribution of the spectroscopic 
sample does not drops at a much brighter limit than that 
of the full photometric sample. These results may be 
somewhat surprising, in the sense that any spectroscopic 
sample tends to be brighter than the photometric sam¬ 
ple in the same survey generally. However in our case, 
we are now focusing only on relatively bright LBGs even 
for the photometric sample, whose magnitude is much 
brighter than the limiting magnitude of the GOSMOS 
survey. This is because we would like to remove most 
contaminants from our sample and select our sample 
whose errors of photometric redshift are as small as possi¬ 
ble. In addition, it is also expected that bright LBGs and 
the quasar sample show a strong correlation because it 
is reported that the LBG clustering becomes strong with 
increasing UV luminosity (e.g., Ouchi et al. 2004), and 
then it is considered that brighter LBGs exist in more 
massive dark matter halos. Therefore bright LBGs are 
useful to investigate the clustering of quasars and galax¬ 
ies. As our measured CCF and AGF will become weaker 
than the real GGF and AGF due to these low- 2 ; galaxies, 
we correct our CCF and AGF by utilizing the contam¬ 
ination rate for color-selected galaxies (see Sections 3.1 
and 3.2). 


Using the samples of qnasars and LBGs described in 
Sections 2.2 and 2.3, we calculate the quasar-LBG two- 
point ang ular CCF, oiqt. (0). at 2 : ^ 4 using the following 
equation (|Croft et al.lll999l : iFrancke et al.ll2008ll : 




(^q-R) 


( 10 ) 


where (DqDi^) and (DqR) are the normalized quasar- 
LBG and quasar-random number of pairs defined as fol¬ 
lows: 


(^Q^l) 


DQDLje) 

VqAl 


( 11 ) 


and. 


{DqR) 


DqR{0) 

VqAr ’ 


( 12 ) 


where DqDi,{9) and DqR{6) are the number of data- 
data and data-random pairs at angular separation 6 ± 
A9, respectively. In the equations (11) and (12), Nq, 
Vl, and Vr are the total number of the quasar, LBG, 
and random sample, respectively. We create the 100,000 
random samples which are avoiding the masked regions 
and we calculate the quasar-LBG CCF. The errors in 
; ^OT,(^) are estima ted by the bootstrap method as follows 
(|Ling et al.lll986fl : 




N(^) - (<^(^))]^ j 


(13) 


where V is the number of bootstrap samples and (w(0)) 
is calculated as follows: 


( 02 ( 0 )) 



2=1 


(14) 


3. CLUSTERING ANALYSIS OF QUASARS AND LBGS 

3.1. Quasar-LBG Two-Point Angular 
Cross-Correlation Function ai 2 ; ~ 4 


We use N = 1, 000, and Figure 5 shows the result of the 
quasar-LBG two-point angular CCF at z ^ 4. Previous 
investigators mentioned that the Poisson error becomes 
increasingly inaccurate at larger scale (e.g., Mountrichas 
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Fig. 3.— Two-color diagram of g' — r' vs. r’ — z\ that we use for 
2 : ~ 4 LEG selection. Orange circles denote the LEG candidates 
with z' < 25.0. Orange circles with a black arrow show the LEG 
candidates which are not detected in ^'-band at 3 sigma limiting 
magnitude, g' = 27.05 (242 objects among 835 objects). For those 
cases, the 3cr lower limit of their g' — r' color is plotted. The 
blue and red lines are the color track of the model star-forming 
galax y (where the instantaneous-burst model of IBruzual &: Chariot! 
120031 with a metallicity of Z = 0.02, an age of 0.025 Gyr and 
a Lyo equivalent width = 21A (| Jones et al.lf20T^ are adopted) 
and elliptical galaxy colo rs (where the stellar population model of 
lEruzual fc Gharlotl 1200"^ with a metallicity of Z = 0.02 and an 
exponential decay time of t= 1 Gyr are adopted, and the ages of 
the model elliptical is 8 G yr). Th e IGM absorption is corrected 
by adopting the model of IMad^ lll995lh for both color tracks. 
TViangles denote the color of the model elliptical and star-forming 
galaxy at 2 : = 1.0, 2.0, 3.0, 3.5, 4.0, and 4.5, resp ectively. Elac k 
asterisks show colors of G, K, and M-type stars JPickle^ 1199811 . 
The black solid line shows our photometric criteria used to select 
LEG candidates at 2 ~ 4. The error bar of the upper right side in 
this Figure denotes Icr error for the g' — r' and r' — z' of our LEG 
candidates. 


et al. 2009). We also confirmed that the errors which are 
calculated by the bootstrap method are about two times 
larger than the Poisson errors. We calculate the CCF 
for both total and spectroscopic samples. We summarize 
wql(^) and errors of a;QL(6*) in Table 2. The observed 
two-point angular CCF is approximated by a power law 
form at large scales as follows: 

ccql( 0 )=A“^( 0 -^-C), ( 15 ) 

where /3 is fixed to be 0.8 (F rancke et al. 2008) an d 
C is the integral constraint (IGroth &: PeeblesI 1197711 . 
To avoid the one-halo term and negative value at the 
second bin, we do not use u}ql{0) on smaller scales 
(angular separation < 80 arcsec) to calculate Here 

we estimate the integral constraint as follows: 


Y.RR{e)e-^- 

ERR{9) 


(16) 



Redshift 

Fig. 4.— Photometric redshift distribution of the 835 color- 
selected galaxies used in this paper. 



Fig. 5.— Spectroscopic redshift distribution of the 67 objects 
among the 835 color-selected galaxies used in this paper. 


We calculate C and a value of C = 0.00601 in this case. 
Using equation ( 15 ) and ( 16 ), we calculate to fit 

the two-point angular CCF. As a result, we obtain 
= 1.77lo gg 3.33l^ gQ for the total and the spectro¬ 
scopic samples, respectively. 

Since we use the photometric sample of LBGs, the de¬ 
rived correlation amplitude is affected by some contam¬ 
ination of galaxies at lower redshifts. Accordingly the 
effect of the contamination on the derived correlation am¬ 
plitude should be corrected. We estimated the contami¬ 
nation rate, /c, utilizing the distribution of spectroscopic 
redshifts in our LEG sample (see Figures 5). In a naive 
estimate, the contamination rate is calculated through 
the following formula using the spectroscopic redshift of 
our color-selected LBGs: /c = {N{zsp < 3 . 1 ) -I- N{zsp > 
4.5))/iVtotai, where the numbers of contaminating objects 
are N{zsp < 3 . 1 ) = 2 and N{zph > 4.5) = 0 , among 
our color-selected and spectroscopically confirmed ob¬ 
jects (iVtotai = 67 ). Therefore the contamination rate 
in our color-selected LBG sample is estimated to be 
fc = 2/67 ~ 0 . 03 . In order to calculate (CCF) ac¬ 
curately, we need to calculate the contamination rate in 
the total sample for low-luminosity quasars. Therefore 
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Fig. 6.— Magnitude distribution of the color-selected galaxies 
used in this paper. Upper and lower panels show the magnitude 
distribution of 835 color-selected galaxies (total) and 67 spectro¬ 
scopically confirmed galaxies at z 4, respectively. 


we calculate the contamination rate for the total sample, 
/cq, as follows: 


/cq = 


^total(-^)[/r^LBG(-^)/*&QSo(-^)]dM 

-Vtotal 


(17) 

where fr, ntotai(M), $lbg(A7), $qso(V/), and Vtotai 
are the fraction which high-redshift galaxies pass quasar 
selections, the magnitude distribution of low-luminosity 
quasars for total sample, the luminosity function of 
LBGs, the luminosity function of quasars, and the to¬ 
tal number of low-luminosity quasars, respectively. Mas¬ 
ters et al. (2012) estimated fr and they found that 
only 2 objects among 386 spectroscopically confirmed 
high-redshift galaxies are pass their quasar selection. 
Therefore /r is estimated to b e ^ 0.005. We use 
‘hLBG('M) which are derived by iBouwens et al.l 1)20071) 
and ^qso{M) which are derived by Masters et al. 
(2012), in this paper. Using above results and Equation 
(17), we calculate /cq and the calculated /cq is ^ 0.05. 
Since the contamination rate for the total sample is so 
low, we do not worry about impact on the CCF due to 
include high-redshift galaxies in the total sample. Even 
if some high-redshift galaxies such as bright LBGs are 
including in the total sample, we consider that the clus¬ 
tering signal will not become weak by the contaminants 
such as bright LBGs because the LBG AGE becomes 
strong with increasing UV luminosity. From the above, 
we do not use the contamination rate for the total sam¬ 
ple though we use the contamination rate for LBGs to 
calculate the (GGF). 

For taking the estimated contamination rate into ac¬ 
count, we calculate A^j (CCF) using the following equa¬ 
tion, 

jCGF 

^-(CCF) = (18) 


We then obtain A^(CCF) =1.82l();gg and for 

the total and the ^ectroscopic samples, respectively. 
The results of and Ai^(CCF) are given in Table 


3. We find that the observed CCF, a;QL(6*) is larger than 
0 at smaller scales. However, there is some possibility 
that this result is only caused by the position of LBGs 
because of the small number of quasars (i.e., this result 
is not caused by the position of quasars). 

In order to confirm that this result (a;QL(0) > 0) is 
caused by the position of quasars and LBGs, we also 
generate 16 random points for quasars (Bq) and calcu¬ 
late the random-LBG CCF 6000 times. Using wrql(0) = 

— 1, we find that a;RQL(6*) 0 at all scale while 

the errors are large at smaller scales (see Figure 7). We 
also calculate the probability of (CCF) for random 
quasars is larger than that of A^, (CCF) for real quasars. 
As a result, this probability is ~ 20%. Since this proba¬ 
bility seems to be high, we treat A(^(CCF) -|-cta,^(gcf) as 
the 86% upper limit. We then calculate the 86% upper 
limits of the spatial correlation length and the bias factor 
for low-luminosity quasars in the same way. 


3.2. LBG Two-Point Angular Auto-Correlation 
Function at z A 

To constrain the quasar triggering mechanism, we have 
to calculate the CCF and the LBG two-point angular 
auto-correlation function (ACF). Since we have calcu¬ 
lated the CCF in Section 3.1, we calculate the LBG 
ACF at z ~ 4 in Section 3.2. In order to calculate 
the LBG two-point an gular ACF at z ^ 4, we use the 
ILandv fc Szalavl (|1993ll estimator: 


wll(^') 


{DD) - 2{DR) -b {RR) 

{m} 


(19) 


where {DD), {DR), and {RR) are the normalized data- 
data, data-random, and random-random number of 
pairs. Those are defined as follows: 


{DD) = 


DlDl{0) 

A'L(fVL-l)/2’ 

DLR{e) 

A'lVr ’ 


{DR) = 


and, 


{RR) = 


RR{e) 


Nn{Nn-l)/2’ 


( 20 ) 

( 21 ) 

( 22 ) 


where DlDl{0), DlR{0), and RR{6) are number of 
data-data, data-random, and random-random pairs at 
angular separation 9 ± A0, respectively. In equations 
(20)-(22), IVl and Ar are the total number of data in the 
LBG and random sample, respectively. Using 100,000 
random samples which are avoiding the masked regions, 
we calculate the two-point ACF of LBGs. The obtained 
results are shown in Figure 7 and Table 4. We find that 
the quasar-LBG two-point CCF is similar with the LBG 
two-point ACF at z ~ 4. The observed two-point an¬ 
gular ACF is also approximated by a power law form at 
large scales as follows: 


a;LL(0)=AACF(0-/5_c), (23) 

where /3 is fixed to be 0.8 (Francke et al. 2008) and C 
(which gives a value of C = 0.00601 in this case) is the 
integral constraint. To avoid the one-halo term, we do 
not use uJi,i^{9) on small scales (angular separation < 40 






























i 



10 100 1000 


Angular Separation (arcsec) 


Fig. 7.— Two-point angular quasar-LBG CCF (red open squares and blue open circles), the random-LBG GCF (black crosses), and the 
LBG AGF (filled black triangles) at 2 ~ 4 in the COSMOS field. Top and bottom panel shows the two-point angular quasar-LBG CCF 
and the LBG AGF in double-logarithmic and semilogarithmic graph, respectively. The data are shown with slight shifts in the horizontal 
direction for clarity. 


TABLE 2 

Summary of the at z A 


6 

(arcsec) 

<0 

(Total) 

DQDCe) 

(Total) 

^QL 

(zsp only) 

(zsp only) 

DQDCd) 

(zsp only) 

14 

0.4970 

0.9997 

2 

1.2178 

1.4805 

2 

46 

-0.1950 

0.2877 

9 

-0.1506 

0.3389 

6 

99 

0.0838 

0.2156 

32 

0.1674 

0.3189 

22 

262 

0.0463 

0.0625 

382 

0.0523 

0.0693 

242 

792 

0.0284 

0.0225 

2803 

0.0347 

0.0282 

1832 


TABLE 3 

Summary of the quasar-lbg ccf 


Nq 

.IVlbg 

QSO Magnitude 

LBG Magnitude 

Redshift 

4 GGF 

-Ajjj 

■A-nj 

ro 








(86% upper limit) 

(/i-i Mpc) 

16“ 

835 

—24.0 < Afi45o < —22.0 

z' < 25.0 

3.1 < z < 4.5 

Q OQ+>i.Ub 

'^•^'^-1.64 

< 6.59 

< 10.72 

25'' 

835 

—24.0 < Mi45o < —22.0 

z’ < 25.0 

3.1 < z < 4.5 

1 rj^j+1.66 

-LA 1 o.86 

i.OZ_Q gg 

< 3.53 

< 7.60 


“Number of the spectroscopically confirmed quasars. 

^Total number of the quasars. 

arcsec) to calculate We calculate the correlation 

amplitude of the LBG AGF, using equations (16) 

and (23). As a result, is S.GSttas- addition, we 

calculate the correlation amplitude for AGF, (AGF) 
as follows: 

jACF 

A. (AGF) = (24) 


the calculated Atj(AGF) is 3.88li;47. These results are 
listed in Table 5. 

3.3. The Spatial Correlation Function 

In order to calculate the bias factor for LBGs and low- 
luminosity quasars, we calculate the spatial correlation 
function for the AGF and GGF. The spatial correlation 
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TABLE 4 

Summary of the ojll at z ~ 4 


e 

(arcsec) 

‘Vll 


Dj^Dj^ie) 

14 

0.8034 

0.5134 

40 

46 

0.1576 

0.1426 

220 

99 

0.0463 

0.0816 

512 

262 

0.0452 

0.0280 

5851 

792 

0.0047 

0.0093 

42091 


function is given as follows: 

^(r) = (r/ro)”'^, (25) 

where vq is the spatial correlation length and 7 = /3+1. 
The spatial correlation func tion for the ACF is calcu lated 
with the following relation (|Totsuii fc Kiharalll969[l . 


A^(ACF) 


I ^{z)Nl{z)E{z)dz 

{c/H,)[! N^{z)dzf 


(26) 


where is the comoving radial distance, Ng{z) is the 
redshift distribution of LBGs, and F{z) is the evolution 
of clustering with redshift, which is assumed to be neg¬ 
ligible and is set equal to 1 in this paper. In equation 
(26), Hj and E{z) are calculated as follows: 


= 


r(i/2) 


r[(7-i)/2] 

r(7/2) 


(27) 


and, 

E{z) = [n^{i + zf+ (28) 


The spatial correlation functio n for the CCF is also c al- 
culated following the relation (|Croom fc Shankslll999ll : 


A fCCF) = '^iz)Ng{z)Nq{z)E{z)dz 

^ {c/Ho)[jNgiz)dzjN^iz)dz] 

(29) 

where Nq{z) is the redshift distributions of quasars. We 
calculate the spatial correlation lengths of the LBG AGF 
and the quasar-LBG GGF using these equations. As a 
result, the spatial correlation length for the LBG AGF is 
6.52t^ gg h~^ Mpc. The 86% upper limits of the spatial 
correlation lengths for the quasar-LBG CCF are 7.60 h~^ 
Mpc and 10.72 h~^ Mpc for the total and the spectro¬ 
scopic sample, respectively. These results are also sum¬ 
marized in Tables 3 and 5. 


3.4. The Bias Factor 

We now investigate the bias factor for low-luminosity 
quasars at 2: ~ 4 in the COSMOS field to study the lumi¬ 
nosity dependence of the quasar clustering and constrain 
the triggering mechanism of the quasar activity. In order 
to investigate the quasar bias factor, we have to estimate 
the galaxy bias at the same redshift. In this paper, we 
estimate the galaxy bias factor at the same redshift us¬ 
ing LBGs. The bias factors for LBGs and quasars are 
defined as follows: 


quasars, and dark matter halos evaluated at 8 h~^ Mpc, 
respectively. Th e correlation f unction of dark matter ha¬ 
los is as follows (lPeebleslll980D : 

Cdm(8,2) = \ (31) 

02 


where J 2 = 72/[(3 — 7)(4 — 7)(6 — 7)2'*'] and al{z) is the 
dark matter density variance in a sphere with a comov¬ 
ing radius of 8 h~^ Mpc. We calculate 0-8(2) using the 
following equation: 


0-8(2) = 08 


Djz) 

74(0)’ 


(32) 


where D{z) is the linear growth factor scaled to unity 
at the present ti me and we calculate it as follows 
(|Mvers et al.ll200^ : 


74(2) 


ffz 1 
90 (1 + 2)' 


(33) 


We calculate by using the fo llowing equation 
(jCarroll et al.lll99^ iMvers et, al.ll2006l l: 


517, 


gz = 




(34) 


where we also c alculate idmz and IIaz as follows 
(|Mvers et al.ll200^ : 


f^mz= (|^)V(l + ^)",f^Az= (35) 

where 77^ is expressed as follows (|Mvers et al.ll2006l) : 


H,=Hopm{l + zf + nA]^^^. (36) 

Using equations (30)-(36), we calculate the correlation 
function of the dark matter halos evaluated at 8 h~^ 
Mpc and the bias factor for LBGs. As a result, the bias 
factor for the LBGs is 4.92^^-29. Next we describe how 
to calculate the bias factor for low-luminosity quasars 
at 2 ^ 4 using the bias factor for LBGs. In order to 
calculate the bias factor for low-luminosity quasars, we 
use the relations as follows (IMountrichas et ^11200911 : 

^LBG = ^8,LBg/'^8,DMj ^QSO^LBG = (37) 

where 6 ql is the bias factor for the quasar-LBG GGF 
and this is calculated as follows: 


^QL — '^8,QSO-LBg/'^8.DM- (38) 

We then obtain the 86% upper limits of 6 ql = 5.65 and 
7.69 for the total and the spectroscopic sample, respec¬ 
tively. Based on these values, we derive the 86% upper 
limits of ^Qgo = 5.63 and 10.50 for the total and spec¬ 
troscopic sample, respectively. The results of the bias 
factors for 6 lbGj ^QLj and bqgo are summarized in Table 
6 . 


&LBg(^:) =-1 


I ^l{8,z) 
^dm(8, 2) 


,bQso{z) =4 


I ^q(8,z) 

Cdm(8, 2) ’ 


(30) 


where ^l( 8,2) = (ro(2)/8 )'>', ^q( 8,2) = (ro(2)/8)'>', and 
Cdm( 8,2) are the spatial correlation functions of LBGs, 


3.5. The Typical Dark Matter Halo Mass 

The inferred bias factor for low-luminosity quasars can 
be used to calculate the typical dark matter halo mass 
(Mdm)- To calculate the typical dark matter halo mass 
in which low-luminosity quasars at 2 ^ 4 reside, we 
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TABLE 5 

Summary of the LBG ACF 


^LBG 

LBG Magnitude 

Redshift 

4ACF' 

■^w 


TO 

(h~^ Mpc) 

835 

z' < 25.0 

3.1 < 2 < 4.5 

3.65+(;^g 

o.oo i 4Y 

^•o^1.96 


TABLE 6 

Summary of the bias factor 


Nq 

Wbg 

QSO Magnitude 

LBG Magnitude 

Redshift 

bQSO 

buBG 

bQL 

16 

835 

—24.0 < A7i45o < —22.0 

z' < 25.0 

3.1 < 2 < 4.5 

< 10.50 


< 7.69 

25 

835 

—24.0 < Afi45o < —22.0 

2 ' < 25.0 

3.1 < 2 < 4.5 

< 5.63 

4.92tl°2l 

< 5.65 
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Fig. 8.— Bias factor for LBGs as a function of UV magnitude. 
Star, square, open circle, triangle, and solid circle points show 
our result, results of lHilde brandt et al.l (12 0091) , lOuchi et al.l 1120041') , 
lAllen et al.l 1120051') . and lOucE^^II.1 1120011) 7 respectively. 


use Equation (8) of iSheth et al.l (1200 It) with parameters 
which are recalibrated bv lxinker et al.l (120051) as follows: 


b{M, z) = l + 


y/aSc l- 


av'^^/a + b\/a{aiy‘^)^^ 


{av'^Y + ^(1 ~ c)(l — c/2)J ’ 


(39) 


where v = , a = 0.707, b = 0.35, c = 0.80, and 

Sc ~ 1.686 is the critical overdensity. a{M) is the linear 
theory rms mass fluctuation on the mass scale M and 
we calculate it using; equations (A8), (A9), and (AlO) of 
Ivan den Boschl (I2002H . The dark matter halo mass can 
be estimated more accurately by utilizing the halo occu¬ 
pation distribution (HOD) models. However the errors 
in the derived quasar bias are still large and therefore we 
only use Equation (39). We then obtain the 86% upper 
limits of log {MDM/{h~^MQ)) = 12.7 and 13.5, for the 
total and the spectroscopic sample, respectively. 


4. DISCUSSION 

We calculated the quasar-LBG two-point CCF and 
the LBG AGE at ^ ^ 4 to investigate the luminos- 



0 2 4 6 


Redshift 


Fig. 9.— Quasar bias factor as a function of redshift. Red, green, 
blue and black squares show the bias factors for our results (total 
and spectroscopic sample), low-luminosity quasars (M 1450 > —24), 
high-luminosity quasars (M 1450 < —24), and high-1-low-luminosity 
quasars, respectively. Black lines show the quasar bias factors as a 
function of redshift with three different models. The model shown 
in the top panel assumes that SMBH growth shuts down after 
quasar phase (efficient feedback). The model shown in the middle 
panel assumes that all 2 > 2 black holes grow with the observed 
quasar luminosity function to the characteristic peak luminosities 
at 2 ~ 2, then shut down (inefficient feedback). The model shown 
in the bottom panel assumes that quasar growth tracks host halo 
growth, even after a quasar episode, until 2 = 2 (maximal growth). 
Dashed, Dash-dotted, and solid lines in each panel show the i'-band 
magnitude ranges at i' < 20.2, i' < 22, and i' < 30, respectively 
(Hopkins et al. 2007). 


ity dependence of quasar clustering and constrain the 
dark matter halo mass in which low-luminosity quasars 
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TABLE 7 

Summary of the reported bias factor for quasars 


Sample 

A^q 

Quasar (AGN) Luminosity 

2 :-Range 

Quasar (AGN) Bias 

Type of CF 

Reference 



(mag or erg s~^) 





COSMOS 

16 

—24.0 < Mi45o < —22.0 

3.1 <2 <4.5 

< 10.50 

CCF 

This work 

COSMOS 

25 

—24.0 < Mi45o < —22.0 

3.1 <2 <4.5 

< 5.63 

CCF 

This work 

GOODS'^ 

25 

—30 < Muv < —25 

1.6 < 2 < 3.7 

3.9 ±3.0'' 

CCF 

Afielberp-er Steide^^fliJi^ 

GOODS 

54 

-25 < Muv < -20 

1.6 < 2 < 3.7 

4.7 ±1.7'' 

CCF 

Afielberper Steide^^fliJ^^ 

2QZ‘= 

18,066 

18.25 <bj < 20.85 

0.3 < 2 < 2.2 

2.02 ±0.07 

ACF 

Croom et al. (2005) 

2QZ 

~ 14,000'' 

Mbj < -22.5^ 

0.80 < 2 < 1.06 

1 CY+U..JU 

•‘•■^'-0.37 

ACF 

Porciani & Norbere (2006) 

2QZ 

~ 14,000'' 

Mbj < -22.5^ 

1.06 < 2 < 1.30 

-•-•'O_o43 

ACF 

Porria.ni Nnrberp- (222^) 

2QZ 

~ 14,000'' 

Mfo^ < -22.5 ^ 

1.30 < 2 < 1.51 

2.33t»;| 

ACF 

Porciani Norbere (2006) 

2QZ 

~ 14,000'' 

Mbj < -22.5^ 

1.51 < 2 < 1.70 

ACF 

Porciani Norbere (2006) 

2QZ 

~ 14,000'' 

Mbj < -22.5 ^ 

1.70 < 2 < 1.89 

O.U2_q 53 

ACF 

Porciani Norbere (2006) 

2QZ 

~ 14,000'' 

Mbj < -22.5^ 

1.89 < 2 < 2.10 

4 -,0+0.49 
+--1-O_0.55 

ACF 

Porciani Nnrberp (2^2^) 

SDSS DRl 

100, 563'' 

14.5 <g < 21.0 

0.4<2<1.0 

1.34 ±0.56 

ACF 


SDSS DRl 

100, 563'' 

14.5 <g< 21.0 

1.0 < 2 < 1.4 

2.20 ±0.26 

ACF 

^^er^^^222£) 

SDSS DRl 

100, 563'' 

14.5 <g< 21.0 

1.4<2<1.7 

2.58 ±0.35 

ACF 


SDSS DRl 

100, 563'' 

14.5 <g< 21.0 

1.7<2<2.1 

2.42 ±0.39 

ACF 

^^er^^^222^) 

SDSS DRl 

100, 563'' 

14.5 <g< 21.0 

2.1 < 2 < 2.8 

3.12 ±0.80 

ACF 

Mver^^L (2006) 

2QZ + 2SLAQt 

22,4169 

18.25 <bj < 20.85'' 

0.3 < z < 2.9 

1.5 ±0.2 

ACF 

da Anpjela et al. (2008) 


6,374* 

20.50 <g < 21.85^ 





MUSYC ECDF-S*^ 

58 

R<25.5 

2.7<2<3.8 

5.5 ±2.0 

CCF 

Francke et al. (2008) 

AEIGS' 

113 

Mb = —20.98 (median) 

0.7<2<1.4 

1.48 ±0.12 

CCF 


XMM-COSMOS"* 

538 

i' < 23.0 

0.2 < 2 < 3.0 

2.0 ±0.2 

ACF 

Gilli et al. (2009) 

2SLAQ 

503 

18.00 <g< 21.85 

0.35 < 2 < 0.75 

1.84 ±0.3 

CCF 

Moimtrichas et al. (2009) 

SDSS DR5 

2,476 

i < 19.1 

0.25 < 2 < 0.60 

1.09 ±0.15 

CCF 

Padmanabhan et al. (2009) 

SDSS DR5 

30,239 

Mi < -22 

0.3 < 2 < 2.2 

2.06 ±0.03 

ACF 


SDSS DR5 

7,902 

Mi < -22 

0.1 < 2 < 0.8 

1.32±0.17 

ACF 

Shen et al. (2009) 

SDSS DR5 

9,975 

Mi < -22 

0.8 < 2 < 1.4 

2.31 ±0.22 

ACF 

Shen et al. (2009) 

SDSS DR5 

11,304 

Mi < -22 

1.4<2<2.0 

2.96 ±0.26 

ACF 

Shen et al. (2009) 

SDSS DR5 

3,828 

Mi < -22 

2.0 < 2 < 2.5 

4.69 ±0.70 

ACF 

Shen et al. (2009) 

SDSS DR5 

2,693 

Mi < -22 

2.9 < 2 < 3.5 

7.76 ± 1.44 

ACF 

Shen et al. (2009) 

SDSS DR5 

1,788 

Mi < -22 

3.5 < 2 < 5.0 

12.96 ±2.09 

ACF 


SDSS DR5 

1,788 

Mi < -22 

3.5 < 2 < 5.0 

9.85 ±2.27 

ACF 


SDSS DR7 

37,290 

i < 19.1 

0.8 < 2 < 2.2 

1.50 ±0.37 

ACF 

Tva.shrhenkn et a, ^^2^^) 

Bootes^ 

445 

log Lbol = 45.86 (median) 

0.7<2<1.8 

2.17 ±0.55 

CCF 


Bootes 

445 

log Lbol = 45.86 (median) 

0.7<2<1.8 

2.50 ±0.65 

ACF 

Jbcko>^f^1^2^^) 

RASS° 

1,552 

43.7 < log Lx (0.1-2.4keV) <44.7 

0.16 < 2 < 0.36 

1.30 ±0.09 

CCF 

Miyaji et al. (20M) 

RASS 

629 

43.05 < log Lx (0.1-2.4keV) < 44.12 

0.07<2<0.16 

1 IQ+O.OS 

ACF 

Krumpe et al. (2012) 

RASS 

1,552 

43.69 < log Lx (0.1-2.4keV) < 44.68 

0.16 < 2 < 0.36 

LOetS;™ 

ACF 

Krumpe et al. (2012) 

RASS 

876 

44.25 < log Lx (0.1-2.4keV) < 45.04 

0.36 < 2 < 0.50 

0.96t“;^5= 

•j qq+012 

ACF 

Krumpe et al. (2012) 

XMM/SDSS 

297 

41.0 < log Lx (2-lOkeV) < 44.0 

0.03 < 2 < 0.2 

CCF 

Moiintrirhas ^ Genrpa.ka,kis ('2(112) 

ManyP 

1,466 

log Lx (0.5-8keV) > 41.0 

0 < 2 < 3.0 

2.26 ±0.16 

ACF 


SDSS DR7 

8,198 

-28.693 < Mi> < -22.576 

0.3 < 2 < 0.9 

1.38 ±0.10 

CCF 

Shen et al. r2013) 


“The Great Observatories Origins Deep Survey (GOODS: IDickinson et ani2003l ). 

*’The bias factors for quasars which ar e derived by Franc ke et al. (2008). 

“The 2dF QSO Redshift Survey (2QZ: [Bovle et al]l20n(1'l . 

'^The total numbers of quasars at 0.80 < 2 < 2.10 llPorcian i &: Norberg|l200 6l'). 

“The absolute magnitude range of quasars at 0.80 < z < 2.10 IIPor'cTaniTr Norbergll200(il L 
iThe 2dF-SDSS LRG and QSO (2SLAQ: IGannon et al.ll200g : [Groom et al.ll2009bl ~). 

®The number of 2QZ quasar sample. 

^The magnitude range of 2QZ quasar sample. 

“The number of 2SLAQ quasar sample. 

■^The magnitude range of 2SLAQ quasar sample. 

^The Multiwavelength Survey by Yale-Ghile (MUSYG; IGawiser et al.||2006l) Extended Chandra Deep Field-South (ECDF-S; ILehmer et al.| 
l 20 ^:rVirani et al.|[ 20 M . 

* All-Wavelength Extended Groth Strip International Survey (AEGIS: IDavis et al.ll2 007ll. 

““XMM-Newton wide-field survey in the COSMOS fiel d (XM M-COSMOS: [Hasmmr et alj|2007l : ICappelluti et alj|2007ll . 

“Bootes multiwavelength survey (Bodtes^JHidgn^tj^ 20071 1. 

“The ROSAT All Sky Survey (RASS: l\^es et alJIlOOflrr 

PChandra Deep Field (CDF)-North, CDF-South, AEIGS, COSMOS, and Extend CDF-South. 


at 2; ^ 4 reside. Since the LBG A CF at z ^ 4 has 
been studie d by many investig a tors llOuchi et alJ 


2004 [20^ lAllen et al.l [200^ iKashikawa et al.l 


20011 


200 


Lee et al.l 120001: Hildebr andt etal.1 120091 : iSavov et alJ 
201 ll : iBarone-Nugent et al.ll2014D . we compare the pre¬ 
vious results of 6lbg with the derived 5lbg in this work 
(Figure 8). Our inferred bias factor for LBGs (6lbg is 
4.92^^ 21 at Muv < —20.9) is not inconsistent with the 
luminosity dependence of the galaxies bias factor which 
are reported from previous studies. We also found that 
the quasar-LBG GGF shows similar clustering with the 
LBG ACF while the errors of our results are large, due 
to the low numbers of low-luminosity quasars. 

To constrain the triggering mechanism of the activity 


in low-luminosity quasars at z ^ 4, we calculated the 
bias factor for low-luminosity quasars in the COSMOS 
field, which are fainter than the characteristic absolute 
magnitude^^. The 86% upper limits of bias factors of 
z ~ 4 low-luminosity quasars are 5.63 and 10.50 for the 
total and the spectroscopic samples, respectively. The 
inferred bias factor for the total sample is smaller than 
that for the spectroscopic sample, while both results are 
just the 86% upper limits. In order to clarify the accurate 
bias factor for low-luminosity quasars, we need larger 

The characteristic absolute magnitude is the absolute magni¬ 
tude where the QLF changes its slope from steep at the brighter 
side to shallow at the fainter side, that is seen typically at M 1450 ~ 
—24 (see e.g., Ikeda et al. 2011). 
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samples of low-luminosity quasars. 

The bias factors for low-luminosity quasars derived in 
this work and results for different redshift and/or lu¬ 
minosity given in previous studies (Adelberger & Stei- 
del 2005; Groom et al. 2005; Myers et al. 2006; Por- 
ciani & Norberg 2006; Ida Angela et al.l 120081 : Francke 
et al. 2008; Gilli et al. 2009; S hen et al. 2009; 
Ross et al. 2009; Coil et al. 2 009; iMountrichas et akl 
120091 : [Padmanabhan et al.l 120091 : Ivachenko et al. 2010; 
Hicko x et al. 2011; Miyaji et al. 2011 ; Krumpe et al. 
2012; IMountrichas fc GeorgakakisI 120121 Koutoulidis et 
al. 2013) are showed in Figure 9. We also summarized 
our result and the previous results of the bias factor for 
quasars in Table 7. The bias factors which are plotted 
in Figure 9 and Table 7 are showing results calculated 
by the ACF or GGF. In order to study the redshift de¬ 
pendence of the bias factor for low-luminosity quasars 
from z ^ 3 to z ~ 4, we compare previous results at 
z ~ 3 with our result at z ~ 4. Francke et al. (2008) 
calculated the bias factor for AGNs at z ~ 3 at the UV 
magnitude range between —26 and —20 and the obtained 
bias factor is 5.5 ± 2. Our result for the total sample is 
consistent with that of their study, though we can not 
rule out the possibility of a lower value of the bias factor 
for low-luminosity quasars. We also compare the bias 
factor for luminous quasars which are brighter than the 
characteristic absolute magnitude at a similar redshift 
(Shen et al. 2009). Shen et al. (2009) calculated the 
bias factor for luminous quasars at 3.5 < z < 5.0 and the 
obtained bias factors are 12.96±2.09 (excluding negative 
data points of the correlation function) and 9.85±2.27 
(including negative data points of the correlation func¬ 
tion), respectively. Our result for spectroscopic sample 
is consistent with their study, though we can not rule out 
the possibility of a lower value of the bias factor for low- 
luminosity quasars. Our result for total sample is much 
smaller than that of Shen et al. (2009). The 86% up¬ 
per limits of the inferred bias factors for low-luminosity 
quasars at z ^ 4 correspond to the typical dark mat¬ 
ter halo mass are log {Mdm/Mq)) = 12.7 and 13.5 
for the total and the spectroscopic samples, respectively. 
This result is not inconsistent with the predicted bias for 
quas ars which is estimate d by the major merger models 
(' e.g.. iHopkins et '^1200711 . 

iHopkins et al.l (|2007ll predicted the bias factors as a 
function of redshift using three different models. The 
first one is that SMBH growth shuts down after each 
quasar phase (efficient feedback). This model is assum¬ 
ing that each SMBH only experiences one phase of quasar 
activity and SMBH growth will stop after this one quasar 
phase. The second one is that all SMBHs at z > 2 grow 
with the QLF to the characteristic peak luminosities at 
z ~ 2 and after that SMBH growth shut down (inefficient 
feedback). This model is assuming that z ~ 6 quasars 
grow either continuously or episodically with their host 
systems until z ~ 2. Hence a quasar feedback at z > 2 
is insufficient to shut down a quasar. The last one is 
that SMBH growth with a Eddington rate until z = 2 
(maximal growth). This model is assuming that SMBHs 
will grow at their Eddington rate until z ^ 2. In case 
of efficient feedback and inefficient feedback model, they 
predict that the bias factor for low-luminosity quasars be¬ 
comes weak with decreasing UV luminosity. In contrast. 


Incase of the maximal growth model, this predict that the 
luminosity dependence of the quasar bias is not detected. 
Our result for total sample is lower than that of the max¬ 
imal growth model. However those three models cannot 
be discriminated by our result for spectroscopic sample. 
While we can also constrain the quasar lifetime using 
the quasar bias factor and quasar space density in prin¬ 
ciple, it is currently too challenging to constrain it due 
to the lack of the spectroscopic sample of low-luminosity 
quasars. In order to measure the bias factor for low- 
luminosity quasars with smaller error bars, we need to 
use larger samples of low-luminosity quasars. Further 
observations of low-luminosity quasars in a wider sur¬ 
vey area are crucial to provide firm constraints on differ¬ 
ent scenarios of quasar evolution and elucidate the trig¬ 
gering mechanism of low-luminosity quasars, especially 
at z > 3, with smaller statistical errors. Surveys for 
high-z low-luminosity quasars with the next-generation 
wide-held prime-focus camera for the Subaru Tele scope 
(Hy per Suprime-Cam: iMivazaki etld] I2006LI20I2I1 , Eu¬ 
clid (iLaureiis et al.ll2011h. an d the Large Synoptic Survey 
Telescope (jivezic et al.ll20d^ will address these issues in 
the very near future. 


5. SUMMARY 

We have estimated the quasar-LBG two-point GGF for 
low-luminosity (—24 < M1450 < —22) quasars and LBGs 
at 3.1 < z < 4.5 in the GOSMOS held. Our quasar 
sample consists of 25 quasars with spectroscopic or pho¬ 
tometric redshifts. This sample is referred as the total 
sample. We also use the 16 quasars with spectroscopic 
redshifts (the spectroscopic sample). We use a sample 
of 835 LBGs with z' < 25.0 in the same redshift range. 
We have also estimated the LBG AGF at z ^ 4 for com¬ 
parison with the quasar-LBG GGF at z ~ 4. Our main 
results are summarized below. 

1. The correlation amplitudes of the quasar-LBG 
GGF are 1.82lQgg and 3.43/j^g4 for the total and 
the spectroscopic sample, respectively. The corre¬ 
lation amplitude of the LBG AGF is 3.88/j^'47. 

2. The 86% upper limits of the spatial correlation 
lengths for the quasar-LBG GGF are 7.60 h~^ Mpc 
and 10.72 h~^ Mpc for the total and the spectro¬ 
scopic sample, respectively. The spatial correlation 
length for the LBG AGF is 6.52l^'gg h~^ Mpc. 

3. The 86% upper limits of the bias factors of z ^ 4 
low-luminosity quasars are 5.63 and 10.50 for the 
total and the spectroscopic sample, respectively. 
The bias factor for the LBGs is 4.92l^'2g. 

4. We find that the bias factor for spectroscopic con¬ 
firmed low-luminosity quasars at z ~ 4 is consistent 
with the bias factor for luminous quasars at z ^ 4, 
though we can not rule out the possibility of a lower 
value of the bias factor for low-luminosity quasars. 

5. We also find that the bias factor for low-luminosity 
quasars at z ^ 4 is consistent with that of previ¬ 
ous results at z ^ 3 at similar quasar luminosity, 
though we can not rule out the possibility of a lower 
value of the bias factor for low-luminosity quasars. 
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6. The 86% upper limits of the inferred dark matter 
halo masses are log Mq)) = 12.7 and 

13.5 for the total and the spectroscopic sample, 
respectively. This result is not inconsistent with 
the predicted bias for quasars which is estimated 
by the major merger models. 

More specific constraints on SMBH growth scenarios 
will be obtained through larger samples of low-luminosity 
quasars at high redshifts that will be discovered by forth¬ 
coming wider and deeper quasar surveys. 
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